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The adsorption behavior of ions at liquid-vapor interfaces exhibits several unexpected yet 
generic features. In particular, energy and entropy are both minimum when the solute resides 
near the surface, for a variety of ions in a range of polar solvents, contrary to predictions of 
classical theories. Motivated by this generality, and by the simple physical ingredients impli- 
cated by computational studies, we have examined intcrfacial solvation in highly schematic 
models, which resolve only coarse fluctuations in solvent density and cohesive energy. Here 
we show that even such lattice gas models recapitulate surprising thermodynamic trends 
observed in detailed simulations and experiments. Attention is focused on the case of two 
dimensions, for which approximate energy and entropy profiles can be calculated analyti- 
cally. Simulations and theoretical analysis of the lattice gas highlight the role of capillary 
wave-like fluctuations in mediating adsorption. They further point to ranges of temperature 
and solute-solvent interaction strength where surface propensity is expected to be strongest. 



2 



I. INTRODUCTION 

Experiments and computer simulations have demonstrated that certain small ions preferentially 
adsorb at liquid- vapor interfaces [IHH] , contrary to expectations from classic theories of ion solva- 
tion [BJ. Efforts to explain this behavior have focused primarily on accounting for effects of solute 
polarizability and the thermodynamic cost of excluding volume in bulk solution (UGHE]- For small 
excluded volumes, this cost is primarily entropic and reflects constraints on the available arrange- 
ments of nearby solvent molecules. [11] . From this perspective the adsorption of a non-polarizable 
ion is expected to be opposed energetically (through loss of dielectric polarization energy) and 
favored entropically (through recovered freedom of local solvent arrangements) . Very recent work 
has shown, however, that these thermodynamic driving forces generally follow the opposite trend, 
with energy and entropy both exhibiting minima when an ion resides near the interface [5j [T2l [T3] . 
In these studies we and others also revealed important driving forces governing ion adsorption that 
had been largely ignored in discussions of interfacial solvation. 

The interface between a dense polar liquid and its coexisting vapor is a region of high energy 
density (relative to that of the bulk liquid), due to incomplete or strained coordination of solvent 
molecules. The liquid's boundary can nonetheless be quite soft, with substantial wavelike shape 
fluctuations even on molecular length scales. These generic facts appear from simulations to have 
important consequences for the thermodynamics of ion adsorption. In particular, moving a so- 
lute from bulk solution to the interface effects more than reduced coordination of the ion (which 
is exclusively emphasized in classic theories). An adsorbed solute also occupies volume in the 
high-energy region, effectively returning solvent density to the bulk environment [3 [T^l 115] . The 
consequence of this simple but often overlooked accounting is a substantially favorable energetic 
contribution. For ions of modest size it is sufficient to render adsorption energies negative [3 [T2T - 
115] . Although electrostatic forces clearly figure prominently in this mechanism, their long spatial 
range is not manifested by the dominant energetic changes accompanying adsorption [5j. Instead, 
a spatially local approximation, which expresses the total energy in terms of the populations of 
different solvent regions (bulk, interface, and solute coordination zone), 

Ulocal — ^bulk^bulk Winter f ace^inter f ace ~i~ £ coot dJl coord ; (1) 

was found to be quite accurate for systems with non-polarizable ions [5]. Here, e 7 is the energy per 
molecule of solvent in region 7, and n 7 is the corresponding number of molecules in that region. 
Softness of the interface, also neglected by conventional theories, is implicated by these simu- 
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lations to significantly influence the entropy of adsorption [5j . A solute that substantially attracts 
solvent, such as an ion in water, tends to suppress surface undulations when it resides near the 
interface. The entropic cost of this suppression can be considerable, opposing adsorption by many 
k^T. Harmonic analysis of capillary wave-like fluctuations in the adsorbed and non-adsorbed states 
was found to roughly account for observed entropy profiles [5]. 

The two physical ingredients speculated in Ref. [3] to be key for ion adsorption, namely surface 
shape fluctuations and heterogeneous energy density near the liquid's boundary, are so generic as to 
be described by the very simplest microscopic models of the liquid state. Motivated by this apparent 
generality, we examine in this paper the classic reduced model of liquid- vapor coexistence, suitably 
adapted to address questions of solvation. Specifically, we consider a two-dimensional lattice gas, 
with a solute that occupies a single lattice cell and attracts only its nearest neighbors. Presenting 
results of both theory and simulation, we demonstrate that the unexpected thermodynamic trends 
revealed by recent atomistic simulations can indeed be captured by such a simple model. The very 
limited microscopic detail in the case of a lattice gas allows unambiguous conclusions regarding the 
source of energy and entropy changes as the solute moves through the interface. It allows as well 
crude predictions for the varying propensity ions exhibit for liquid-vapor interfaces. 

The rest of the paper is organized as follows. In Sec. [TTJ, we describe our lattice gas model for 
solvation. We present results of Monte Carlo simulations in Sec. |III| In Sec. |IV| we develop a theory 
for the behavior of this model based on the so-called solid-on-solid simplification and compare its 
predictions with numerical results. We conclude in Sec. [V] with a discussion of specific implications 
for three-dimensional molecular systems. 



We consider a two-dimensional L x L lattice gas system with a single solute, that occupies one 
lattice cell. Cells are indexed by a horizontal coordinate x ranging from 1 to L and a vertical 
coordinate y ranging from — [L — 1) /2 to (L — l)/2 Q The x-coordinate of the solute is constrained 
to x = (L + l)/2. The solute interacts only with the 4 nearest neighbor cells. It additionally 
excludes volume, so that a cell may not be simultaneously occupied by both solvent and solute. 
This system is described by the Hamiltonian 



II. LATTICE GAS MODEL 
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1 Without loss of generality, we assume that L is odd. 
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FIG. 1. Profiles of (a) energy E(h), (b) entropy S(h), and (c) free energy F(h), as a function of the height 
of the solute, h, for various values of ei at T = eo/3kg. The height of h — coincides with the GDS. These 
profiles were obtained from Monte Carlo simulations of the lattice gas. 



where 

H = -e ^2 rijUj - (j, y~] rij (3) 

i,j&nn i 

describes the pure solvent. The occupation variables n and nf are binary: rtj = 1 if the i th 
lattice site is occupied by a solvent and zero otherwise, nf = 1 if the i th lattice site is occupied 
by a solute and zero otherwise. Y^ijenn denotes a sum over all nearest neighbor lattice pairs i, 
j. The coupling constant e$ > describes attraction between neighboring solvent cells; e\ is the 
solvent-solute coupling constant. 

For both theory and simulation we work with a grand canonical ensemble, in which the number 
of solvent cells can fluctuate subject to chemical potential /i. (The number of solute cells is 
fixed.) Conditions of liquid-vapor coexistence are imposed by setting [i = — 2eo, and applying 
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boundary conditions that favor liquid in y < and vapor in y > 0. Specifically, lattice cells with 
y = — [L— 1)/2 are constrained to be occupied by solvent; those with y = (L — l)/2 are constrained 
to be unoccupied (i.e., vapor). 

The mean interface position, or the Gibbs Dividing Surface (GDS), is fixed at y = by con- 
straining as well the state of lattice cells with x = 1 and x = L; those with y < are occupied, 
and those with y > are unoccupied. 

We will examine values of e\ between and 1. By symmetry (see Appendix A), a solute 
with ei = eo/2 is equally likely to reside in a cell well within the liquid phase as in a cell well 
within the vapor phase. Judging by the free energy of transfer from liquid to vapor, one would 
therefore consider solutes with e\ < eo/2 to be solvophobic, and those with e\ > eo/2 to be 
solvophilic (although the latter may attract solvent significantly less strongly than adjacent solvent 
cells attract each other). 

Using theory and simulation, we have determined the free energy F, average energy E, and the 
entropy S as functions of the height h of the solute (i.e., its y-coordinate). 

III. MONTE CARLO SIMULATIONS 

We present results of Monte Carlo simulations at temperature T = co/3/cb (well below the 
critical temperature T c ~ 6o/1.764A;b) for a lattice with linear dimension L = 30. Trial moves 
included (i) changing the state of a cell from solvent to vapor, (ii) changing a cell from vapor to 
solvent, and (iii) translating the solute along the y axis from a lattice cell s to a cell s' . In the 
case of solute displacements, two different trial moves were used, distinguished by their treatment 
of solvent occupation in the cell s vacated by the solute. In one move, cell s adopts the previous 
state of s' , i.e., n* rm ' = n s i. In another move, cell s adopts the complementary state of s', i.e., 
n triai _ \— ngl t i n oil cases detailed balance was ensured through appropriate Metropolis acceptance 
criteria, preserving the Boltzmann distribution p({nj,n|}) oc e~^ H , where j3 = is inverse 

temperature. 

We obtained estimates of F(h) by computing the equilibrium probability p(h) oc exp[— f3F(h)] 
that the solute resides at height h. The average energy E(h) was estimated directly from simula- 
tions, while estimates of the entropy S(h) were calculated from TS(h) = E{h) — F(h). 

The profiles of F(h), E(h), and S(h) obtained from simulations with ei/eo = {0.3, 0.4, 0.5, 0.6, 0.7, 0.8}, 
are plotted in Figs. [TJ These results echo several nontrivial thermodynamic trends observed in 
molecular simulations. In particular, free energy profiles for e± = {0.4eo, 0.5eo, 0.6eo, 0.7eo} exhibit 
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minima near h = 0. These solutes adsorb to the interface, i.e., their density at h = is higher 
than in either bulk phase. Entropy minima at the same values of ei indicate that this adsorption 
reduces the diversity of accessible solvent configurations. The driving force for surface propensity 
is instead energetic, as evidenced by strong minima in E{h) near h = 0. Interestingly, adsorption 
is strongest for the case ei = eo/2, at the crossover between solvophobic and solvophilic regimes. 

These trends were also observed in simulations with smaller and larger system sizes (L = 
{20,30,40,50,60,70}), and for various temperatures in the range eo/5/ce < T < eQ/2ks- In all 
cases, profiles of free energy, energy, and entropy exhibit minima when the solute is at the interface, 
and maximum adsorption occurs at e\ = eo/2. 
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FIG. 2. Illustration of three configurational subensembles in the RSOS model. The fully solvated in state 
(a) features a liquid vapor interface entirely above the solute. In the pinned state subensemble (b), this 
interface directly contacts the solute. In the unsolvated state out state, (c), the interface lies below the 
solute. 



IV. THEORY 



The 2-d lattice gas model of solvation described by Eq. [2j within approximations that are 
well justified for T <C T c , is sufficiently simple to permit an analytical solution. This section 
presents our mathematical analysis, which draws heavily on results from theoretical studies of 
surface roughening in Ising-like models. [16H19] Our key assumptions are that the liquid- vapor 
interface is well defined and that the solute does not significantly alter fluctuations within either 
phase. More specifically, we restrict attention to the so-called "solid-on-solid" (SOS) limit [20| [2Tj . 
in which: (1) all cells in the liquid phase are occupied by solvent (except the cell occupied by 
the solute when it is present), (2) all cells in the vapor phase are unoccupied (except the solute 
cell), and (3) the boundary between liquid and vapor is a single- valued function of x. With this 
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neglect of bubbles, islands, and overhangs, a configuration can be specified through the heights of 
L columns representing the instantaneous domain of the liquid phase. 

Under conditions of coexistence, the Hamiltonian describing the solvent is then given by [201 122] 

^o = | £ \ht-hjl (4) 

i,j£nn 

where hi denotes the height of the liquid in the i th column. Boundary conditions described in 
Sec. |n] fix the heights of the first and last columns, hi = and hi = 0. They also bound the 
range of column height fluctuations, — L/2 < hi < L/2. As an added simplification, we work with 
the Restricted SOS (RSOS) model [T7] where the height difference between neighboring columns is 
constrained to be at most 1, \hi — hi+i\ < 1. A few allowed configurations of the RSOS model are 
illustrated in Fig[2] Previous work has shown that fluctuations of the neat liquid- vapor interface of a 
two-dimensional lattice gas are well described by the RSOS model for temperatures T <C T c (T7][T9]. 
As in that work, the restriction on height fluctuations is necessary to derive simple closed form 
expressions for thermodynamic quantities of interest. 

The liquid-vapor interface of the RSOS model supports shape fluctuations akin to capillary 
waves. It has been shown in particular that long- wavelength Fourier modes are Gaussian distributed 
to a good approximation. For a 2-d RSOS interface of length M, whose boundaries are fixed at 
hi = p and hu = Q, the partition function is given by [TTJ [191 123] 



JlaM 



-(P ~ g) 
4aM 



(5) 



in the limit Ma >> 1, where a = exp(— /3eo/2). This result is strictly valid only in the low tempera- 
ture limit, where terms of order 0(a 2 ) can be neglected. Empirically, however, the thermodynamic 
properties we compute using Eq. [5] agree well with those obtained from lattice gas simulations even 
at moderate temperatures. 

To construct a theory for interfacial solvation from Eq. [5j we first partition the ensemble of 
solvent fluctuations consistent with a given solute height h. We classify configurations in which 
the interface lies entirely above the solute as belonging to a completely solvated class, denoted in. 
Immersion in this sense requires h s > h, where s is the solute's x-coordinate. Configurations in 
which the solute makes no contact with the liquid phase, h s < h — 1, we classify as out. In the 
remaining configurations, which we classify as pinned, the solute resides at the interface, h s = h—1, 
in effect pinning local height fluctuations. 

Defining Zi n (h), Z out (h), and Z p i nne d(h) as the partition functions for these three subensembles, 
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we can reconstruct the total partition function for a given solute height simply as 

Z(h) = Z in (h) + Z out (h) + Zpi nned {h) . (6) 

Each of the three configurational classes corresponds to a constrained ensemble of pure liquid- 
vapor configurations, reweighted by solvent-solute interactions. We will calculate their partition 
functions in turn. 

The out subensemble is simplest, since solvent and solute do not interact. The partition function 
for this case is given by 

Zout(h)= E---E e ~ Pn °= E Z 1 ' s (0,h s )Z s ' L (h s ,0), (7) 

/i 1= ,h L =0 ,h s <h-l h a <h-l 

where the notation J2" ' J2h!=o h L =o h s <R-i i m phes a sum over all configurations of the RSOS 
system subject to the constraints {hi = , hi, = , h s < h— 1}. To obtain a closed-form expression, 
we convert sums to integrals in the limit L » 1, and set the solute's horizontal position at 
s = (L + l)/2, yielding, 



2aL rh-1 2aL rh-1 

Zout(h) - tt~7 / dh s e~ h ^ « J_ / dhs e-^L 



27roL ./-(l-i)/2 " 2iraL j_ OQ 

m' ( /( ,' ) ' 

aL 




(8) 



Computing the energy Hi n of a fully solvated configuration (belonging to the in state) requires 
accounting for direct interactions between solvent and solute, as well as counting broken solvent- 
solvent bonds and the thermodynamic consequences of removing a solvent cell to accommodate 
the solute. The four solvent-solute bonds contribute an energy — 4e\. Loss of 4 solvent-solvent 
interactions is partially offset by the reversible work [i = — 2eo of returning one solvent cell to the 
bath, yielding a total energy 

H in = H - 4(ei - e /2) . (9) 
The corresponding partition function can therefore be written as 

Z m {h)= e-^e^-^ (10) 

/i 1= 0,/ii,=0,/i s >/i 

Again, using Eq. [5] in the limit L >> 1, and rewriting the sums as integrals, we obtain 



J2aL r (L-l)/2 J2aL r°° 

V ; 27TaLj- h 2itaLh 



4 \VaL 



(11) 
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Statistics of the pinned state are complicated by variation in the solvent-solute interaction 
energy. Within the RSOS model, allowed heights of columns s — 1 and s + 1 are h, h — 1, and 
h — 2. Only in the first case does the solute interact with an adjacent liquid column. Since it is 
not necessary to break solvent-solvent bonds or remove solvent cells in this case, the corresponding 
energy function can be written as 

Hpi nne d = Ho + AH(h s ^i,h s+ i;h), (12) 

where 

AH(h s -i,h s+1 ; h) = -ei(l + S hs _ uh + S hs+1 ~ h ). (13) 

Summation over the column heights fi2, hs, ■ ■ ■ , fo s -2 and h s+ 2, h s +3, . . . , may now be per- 
formed using Eq. [5j 

h h 

Z pi nned= £ ^-^(0, h s ^) exp[- (3 AH (h s ^, h s+1 ;h)]x 

h s -!=h-2 h s+ i=h-2 (14) 

xexp =^(1^-^-1)1 + 1^-^-1)1) z( s+1 > L \h s+1 ,0). 

The constraint h s = h — 1 defining the pinned state, together with the locality of inter-column 
interactions, decouples fluctuations to the left and right of the solute. As a result, summations 
over h s -i and h s+ \ can be carried out independently; with the choice s = (L + l)/2 they yield 
identical contributions in the L>1 limit. We then obtain 

p 2aL 

Z pmned (h) « 6 2 — e-^-D l aL , (15) 

where 

b _ e -2a r e /3(3e 1 /2-e /2) e (-2ft+l)/(2aL) + ^(ei/2) + ^((e, -e )/2) e (2h-3)/(2aL)l ^ ^gj 

Using Eqs. 6|8|ll|l5 we can readily obtain profiles of free energy F(h) = — h\Z{h), energy 



E(h) = —(dlnZ(h)/df3)L, and entropy S(h) = (E — F)/T. These analytical results are plotted in 
Fig.[3]for a lattice of size L = 30, at temperature T = eo/3.6/cB> and for various values of ei/eo- They 
agree very well with simulation results. In particular, the depth of minima at h = are accurately 
predicted. In Fig. |l]we focus on these adsorption properties, Ag a d s = g(0) — g{—(L — l)/2) for 
thermodynamic property g, for the specific case e\ = eo/2^} Theoretical predictions are plotted 
as functions of temperature at L = 30 alongside results of Monte Carlo simulations. 



2 Recall that the bulk values of F, E, S are equal in the liquid and vapor states when ei = eo/2. Hence, the results 
would not have changed if the bulk vapor was chosen as the reference state. 
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h h 
(c) ei = 7e /12 (d) t x = 86 /12 

FIG. 3. Profiles of the energy, entropy, and the free energy as a function of h for various values of ei at 
T = eo/3.6fcs. The blue lines were calculated from Eq. [6j while the dashed black lines with squares were 
obtained from simulations. 



Agreement is very good for temperatures in the range eo/^ks <T < eQ/3ks for L = 30. 

The failure of our theory at high temperature stems from neglect of contributions to Z from 
terms of order 0(a 2 ) and higher. At high temperature ruggedness of the interface can also render 
the RSOS model a poor proxy for the SOS model. At higher temperatures still, neglect of bulk 
inhomogeneities in the SOS model is itself not justified. At low temperatures a different prerequisite 
for the validity of Eq. [5] is violated, namely, aL is no longer a large number. For the RSOS fluid 
with L = 30 columns, aL ~ 1 when T ~ 0.147eo/fe^. For temperatures around or lower than this 
value, height fluctuations in the RSOS fluid are no longer Gaussian to a good approximation. 

Success of the SOS approximation in the temperature range e^/Sks < T < eo/3/ce reveals 
clearly the source of entropy variations observed in our lattice gas simulations. Absent density 
fluctuations in the liquid or vapor phase (which are neglected in the SOS description), fluctuations 
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FIG. 4. Free energy, energy, entropy of adsorption as a function of T for ei = eo/2, with L = 30. The red 
line is calculated from theory, while the black line with circles is obtained from simulations. 



in interfacial shape provide the sole source of entropy in these schematic models. The minimum in 
S(h) near h = thus directly indicates suppression of capillary fluctuations when the solute resides 
at the liquid-vapor interface. Similarly, local maxima of S(h) just above and below h = signify 
enhancement of surface undulations when the solute pulls the interface a few lattice spacings away 
from its natural equilibrium position. 

Interfacial pinning is strongest for a solute with coupling constant e\ ~ eo/2. In simulations 
this conclusion is suggested by the depth of entropy minima at h = 0, which are most pronounced 
at intermediate coupling. It can be drawn more directly in our RSOS theory from the statistical 
weight P p inned(h) = Z p i nne d{h) / Z (h) of the pinned subensemble. For a given value of e\, P p i n ned(h) 



is always greatest at h = 0, where pinning does not require deforming the interface. (See Fig. 5(a) ) 
Ppinned(fy should clearly increase with e\ in the strongly solvophobic regime, where weak attraction 
to solvent is insufficient to offset the entropic price of suppressing capillary-like waves. It is also 
clear that P p i nn ed(0) should decrease with e\ in the strongly solvophilic regime, where extremely 
favorable solvation encourages the fluctuating interface to deform above a solute at h = (i.e., to 



remain in the in state) . Fig |5(b) bears out the resulting expectation that pinning is most effective 
at intermediate ei, very close to e\ = eo/2. 

Subensemble weights P a {h) = Z a (h)/Z(h), where a = {in, out, pinned} , also shed light on the 
nature of entropy maxima at values of h a few lattice spacings above and below the GDS. In the 
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RSOS model the total entropy can be decomposed by subensemble: 

S(h)=J2 S *(h) p *(h)-k B J2P«W^P a (h), (17) 

a a 

where S a = (E a — F a )/T characterizes fluctuations within each subensemble. The final term in 



Eq. 17 is an entropy of mixing associated with fluctuations between in, out, and pinned states. 



Plots of P a (h) in Fig. 5(a) show that in and pinned states become about equally probable near 
h = —2, contributing a mixing entropy of roughly ZcBln2. This gain is partially offset by the 
intrinsically low entropy of the pinned state, and by decreasing Sin as the requirement h s > h 
becomes a nonnegligible constraint on interfacial fluctuations in the solvated state. The residual 
increase in entropy is a signature of spontaneous switching from the pinned state to an unpinned 
state, and vice versa. Analogous switching between pinned and unsolvated states is responsible for 
the entropy maximum at h = +4. 
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FIG. 5. (a) Probabilities of the pinned, in, and out states as a function of h for T — eo/3.6fcs, ei/eo = 1/2 
and (b) probability of the pinned state at h — as a function of t\ for T = eo/3-6fcs. When the solute is 
near the GDS at h = 0, most of the configurations are in the pinned state. 



Energy can also be decomposed by subensemble: 



(18) 
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where E a (h) is the average energy in subensemble a when the solute is held at h. Eqs |4j [9 
allow rough estimates of E a {h) and thus simple predictions for energies of adsorption. For h S> 
the unsolvated state dominates (P ou t = 1), giving E(oo) w E out ~ Eq, where Eq is the average 
energy of the solute-free system. For h <C the solvated state dominates (Pj n = 1), giving 
E(—cc) m Ei n sa Eq — 4(ei — eo/2) from Eq. M To the extent that pinning dominates at h = 0, 
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FIG. 6. Adsorption propensity denned as z(-(l-i)/2)^+,z((.l-i)/2) as a function of ei for T = eo/3-6/cs. These 
estimates were calculated using Eq. [6] 



we have E(0) « E pinne< }{Q) ~ So + AH(h a -i, h s+ i; h). At temperatures low enough to safely set 
fo s _ 1 = /is+i = /i for the pinned state, we have -E(O) £" — 3ei — eo- Interfacial adsorption of 
the solute from the liquid phase therefore produces an energy change E(0) — E(—oo) ei — eo, 
accounting roughly for replacing a solvent-solute bond with a solvent-solvent bond. For a modestly 
solvophilic solute (eo/2 < ei < 1), adsorption is energetically favorable, and the energetic driving 
force strengthens with decreasing e\. This force is stronger still in the solvophobic regime (ei < 
eo/2), but these solutes reside primarily in the vapor phase. The strongest relevant energy of 
adsorption from the liquid phase is thus obtained for e\ ~ eo/2, for which AE a d s ~ —eo/2. 

These crude predictions are consistent at low temperatures (below T = 0.15eo/&B for L = 30) 
with results from our simulations and detailed theoretical analysis (see Fig. [4b. At temperatures 
above T = 0.15 eo/fcs for a lattice of size L = 30, interfacial fluctuations modify AE a d s significantly. 
Reduced coordination at a rough interface enhances the energetic driving force for adsorption, 
effectively increasing the number of solvent-solute and solvent-solvent bonds exchanged when the 
solute moves to the surface (see Fig. [4]). At the same time, however, the opposing entropic forces 
discussed above become more pronounced as temperature increases. It is not obvious from our 
simple estimates how this competition plays out, though the sign and magnitude of AS a ds ensures 
that net surface propensity decreases smoothly with increasing temperature (see Fig. [4]) . Regardless 
of temperature, adsorption is always strongest when t\ = eo/2 (see Fig. [6]). 
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V. CONCLUSIONS 

Among models that resolve microscopic fluctuations of the liquid-vapor interface, the two- 
dimensional lattice gas we have studied may be the simplest from which mechanisms of solute 
surface propensity might be learned. Any insight it offers into the behavior of real molecular 
systems is no doubt schematic. But the fact that it captures thermodynamic trends which elude 
celebrated theories of solvation encourages a detailed view of its conclusions in the context of more 
realistic molecular models. 

Any such comparison requires assigning interaction parameters eo and ei to systems with much 
more complicated energetics. We do so through the local approximation of Eq. [TJ which for 
simulations accurately describes the potential energy of a non-polarizable ion in water at a level 
of detail akin to that of a lattice gas. For a three-dimensional lattice gas model of bulk solvent, 
the total energy could be calculated by attributing to each liquid cell half of the interactions with 
its 6 neighboring liquid cells, so that ebulk = — 3eo (neglecting solvent density inhomogeneities) . 
Similarly, since solvent cells at an ideally flat liquid-vapor interface engage in one fewer interaction 
than in bulk, we estimate ^interface. ~ — (5/2)eo- From the value of ^interface — tbulk = 6.52 kJ/mol 
computed from atomistic simulations of water [5], we thus infer a coarse-grained solvent-solvent 
attraction strength eo ~ 13.04 kJ/mol. Within the local approximation, solvent cells adjacent to 
the solute would be assigned half of their solvent-solvent interactions as well as a single solute bond. 
For a solute well within bulk solution, we then have e coor d = —e\ — (5/2)eo in the absence of solvent 
density inhomogeneities (i.e., within the SOS approximation). From the value e coor( i — e^uik = 
—AAA kJ/mol determined for a fractionally charged non-polarizable Iodine ion in water, I -0 ' 8 , by 
simulation [5J we finally obtain ei = 10.96 kJ/mol, or ei/eo ~ 0.84. 

From a lattice gas perspective, the simulated solute I -0 - 8 (aq) therefore falls within the hy- 
drophilic regime ei/eo > 1/2. Further reducing the magnitude of the solute's charge should cer- 
tainly weaken solute-solvent attractions. Decreasing e\ by a small amount in this way should bring 
the solute closer to conditions of maximum adsorption, ei = eo/2, where the scales of adsorption 
energy and entropy are both largest. Based on results we have presented, we thus expect that 
changing the ion's charge from — 0.8e to — 0.75e would deepen the interfacial minima of energy and 
entropy profiles. Using the same molecular simulation techniques described in Ref. [5], we find this 
prediction to indeed hold true. 

In drawing this comparison we have assumed that basic interfacial solvation properties of the 
lattice gas are insensitive to dimensionality d. Although scaling behavior of a lattice gas near the 
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critical point is famously sensitive to d, there is good reason to expect that the essence of our 
conclusions regarding solvation will be unchanged when moving from d = 2 to d = 3. The liquid- 
vapor interface of a three-dimsional lattice gas similarly supports capillary wave-like fluctuations 
for temperatures in the range Tr < T < T c , where Tr denotes the so-called roughening transition 
temperature [24J. Furthermore, these surface undulations are Gaussian distributed to a good 
approximation at long wavelengths |16| I25| . just as in Eq. [5] Partition functions of the in, out, 
and pinned states should therefore have similar forms in d = 3. In qualitative terms, only the scale 
of height fluctuations governing Z{ n , Z out , and Z p i nnec [ will change, growing with system size only 
as vhiX rather than y/L [22] • Preliminary simulation results for the d = 3 case bear out these 
expectations. 

We conclude that lattice gas models of interfacial solvation are a promising starting point for 
scrutinizing essential effects of surface roughness and energetic heterogeneity on solute adsorption. 
For the especially interesting case of charged solutes in highly polar solvents, a more realistic 
treatment of electrostatics is clearly in order. Given the success of lattice models for dielectric 
response [26], this improvement appears quite feasible. An ability to treat solutes of varying 
size and geometry would also greatly aid connections with experiment. These developments are 
underway. 



Appendix A: Appendix 

Here we present arguments underlying the assertion that a solute with coupling constant e\ = 
eo/2 is equally likely to reside in either bulk phase. Consider a pair of configurations {ni,nf\ 
and {n'^nf} related by a global transformation that swaps liquid and vapor phases. Specifically, 
the configuration {n^,n|} is obtained from {nj,n|} by changing the solvent occupation state of 
all lattice cells not occupied by the solute, n\ = (1 — n|)(l — rij). Repeating this transformation, 
Hi = (1 — n|)(l — nfj), we recover the original configuration {rii,nf}. In simple terms, if the solute 
resides in the liquid phase in {m,nf}, then it resides in the vapor phase in {n' t , nf }, and vice versa. 

With an accompanying transformation of the solute-solvent coupling e\, the energies of config- 
urations {ni,nf} and {n^nf} can be directly related. Defining E({ni, nf}; {eo, ei}) as the energy 
of configuration {n,i,nf} with solvent-solvent coupling eo, and solvent-solute coupling e\, we have 
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£({n-,nf};{e ,ei}) = -e ^ n -n^- - // ^ n- - ei ^ nfn^ 

i,j£nn i i,j£nn 

= -eo ^ n^nj - /x ^ nj - (e - ei) ^ nfn, + 4(ei - e /2) (Al) 

= E({nj,n|};{eo,e - ej) + 4(ei - e /2). 

When ei = cq/2, the energy of the two configurations are the same. Thus for every configuration in 
which the solute is in the liquid phase, there exists a configuration with the same statistical weight 
in which the solute is in the vapor phase and vice versa. Hence solutes with ei = eo/2 exhibit equal 
preference for the liquid and vapor phases. 
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